Pairwise summation approximation for Casimir potentials and its limitations 



(N 

O 

(N 

o 

O: 



> 

00 
(N 
00 

d 

(N 



X: 



Annc-Florence BitbolIB Antoinc Canaguier- DurandQ Astrid Lambrecht, and Serge Reynaud 
Laboratoire Kastler Brossel, ENS, UPMC, CNRS, F-75252 Pans, France 
(Dated: November 1, 2012) 

We investigate the error made by the pairwise summation (PWS) approximation in three ge- 
ometries where the exact formula for the Casimir interaction is known: atom-slab, slab-slab and 
sphere-slab configurations. For each case the interactions are calculated analytically by summing 
the van der Waals interactions between the two objects. We show that the PWS result is incorrect 
even for an infinitely thin slab in the atom-slab configuration, because of local field effects, unless 
the material is infinitely dilute. In the experimentally relevant case of dielectric materials, in all 
considered geometries the error made by the PWS approximation is much higher than the well- 
known value obtained for perfect reflectors in the long-range regime. This error is maximized for 
permittivities close to the one of Silicon. 

PACS numbers: 34.35.-|~a, 12.20Ds, 03.70.-|-k 



I. INTRODUCTION 

The Casimir effect is the universal attraction between 
two perfectly reflecting plates in quantum vacuum due 
the zero-point fluctuations of the electromagnetic field as 
first described in the seminal paper by Casimir—. A quan- 
tum mechanical treatment of van der Waals forces includ- 
ing the finite speed of light had already been established 
by Casimir and Polder shortly before by means of quan- 
tum electrodynamics^. It is therefore natural to wonder 
whether Casimir interactions between macroscopic ob- 
jects can be deduced from the van der Waals interactions 
between their constituting atoms. A simple way of mak- 
ing this link would be to sum the van der Waals interac- 
tions between all pairs of constituting atoms, an approxi- 
mation method known as pairwise summation (PWS). It 
is worth noting that attempts at explaining macroscopic 
interactions between colloids by pairwise summing van 
der Waals interactions were already carried out by de 
Boer and Hamaker— ten years before the work of Casimir 
and Polder. 

However, the attractive PWS idea, that would reduce 
the Casimir eff'ect - a physical effect of boundary condi- 
tions on the electromagnetic vacuum field - to a macro- 
scopic resultant of two-body interatomic forces, has not 
been fundamentally successful. Indeed, PWS fails by 
construction to take into account any many-body effect. 
Such effects, for instance the screening of electromag- 
netic fields, are crucial in condensed matter physics, and 
have to be accounted for when computing the Casimir 
effect between macroscopic objects. Another fundamen- 
tal cause of the insufficiency of PWS is the existence of 
three-body van der Waals-type interactions^. For per- 
meable materials even the sign of the interaction can be 
erroneous when calculated within PWS^. Nevertheless, 
given the complexity of performing exact calculations of 
the Casimir effect in different geometries, the PWS ap- 
proximation is still used nowadays since it provides ap- 
proximations to the true Casimir interactions^^—. The re- 
sults obtained by PWS are then often empirically renor- 
malized by a corrective factor computed in a simpler case 



where the exact Casimir interaction is known. 

Limitations of PWS have been given in the past us- 
ing a path integral formulation for fluctuation-induced 
forces to study the orientational dependence of PWS^, 
the force between deformed plates^ or within a pertur- 
bative expansion in the dielectric contrast between ar- 
bitrarily shaped bodiesii. Recently first principle cal- 
culations have obtained the PWS approximation as an 
asymptotic approximation of the Casimir energy in the 
weakly-coupledi^ or diluted limii^i^. In this paper, we 
investigate the error made by PWS using the scattering 
approacbi^"— in the atom-slab, slab-slab, and sphere- 
slab geometries, three fundamental geometries where the 
exact formula for the Casimir interaction is known. We 
investigate which effects cause the failure of the PWS ap- 
proximation, and we study its error as a function of the 
material and the slab thickness. In Sec. |lll the interac- 
tions between an atom and a slab, between two slabs and 
between a slab and a sphere are calculated analytically 
by PWS at any distance between the two interacting ob- 
jects. Then, in Sec. Illli the importance of local field 
effects as a cause of the failure of PWS is emphasized. 
The influence of the material for bulk mirrors and the ef- 
fect of slab thickness in the long-range limit are studied 
in Sec. IIVI We summarize our flndings in Sec. |Vl 



II. PWS CALCULATION OF CASIMIR 
INTERACTIONS 

PWS calculations are usually carried out in the 
short-range and long-range limiting regime o^'^d'^' J^ where 
Casimir-Polder interactions reduce to power laws. For 
the sake of generality, we carry out the PWS calculation 
of the general Casimir-Polder interaction in order to ob- 
tain PWS estimates of the interactions between an atom 
and a slab, between two slabs, and between a sphere and 
a slab, at any distance between them. In these three well- 
known geometries, the PWS results can then be com- 
pared to the exact ones. 
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A. Atom-slab geometry 

We start from the van der Waals formulap^ giving the 
interaction energy between two atoms A and B of polar- 
izabilities aA{i^) and ^^(aj), lying at a distance d: 



he r 



I'^cP ud 



(1) 



where the integral is written over imaginary frequencies, 
w = icu. 

We then consider a slab of thickness ca constituted of 
atoms per unit volume at a distance L along the z 
axis from an isolated atom B (see Fig. [T]). 



following interaction potential within PWS: 
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B. Slab-slab geometry 



In order to compute the PWS estimate of the Casimir 
interaction between two slabs, we simply sum the pre- 
vious result ([3]) over atoms B constituting another slab 
with thickness and number of atoms per unit volume 
in the same way as before (see Fig. [5]). 
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FIG. 2. Slab-slab geometry. 



FIG. 1. Atom-slab geometry. 

We compute the PWS interaction between the slab 
(s) and the atom (a) by summing the van der Waals 
interactions ([T|) between each atom A of the slab (s) and 
atom B: 



C/r_T(i,e^) = 27r< 



L+eA 



dz / drrUa-aid) (2) 
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where d = \/z^ + r^. The integrations over r and z can 
be carried out by parts after inverting them with the 
integration over u, yielding the following result: 



u!Z'iL,eA)^hcn^ [ 
Jo 



du u 
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x[f{2uL)- f{2u{L + eA))] (3) 

r being the incomplete Gamma function F : (a, x) t-j- 

A limiting case is the infinitely thick slab {ca — ^ oo), 
corresponding to a bulk material. We will call this sit- 
uation atom (a)- bulk plate (p) configuration in the re- 
mainder of the paper. For this geometry we obtain the 



This procedure leads to the following PWS estimate of 
the Casimir interaction between two slabs of thicknesses 
eA and eg at a distance L, per unit surface: 
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x[g{2uL)+g{2u{L+eA+eB))-g{2u{L+eA))-g{2u{L+eB)) \ 



with 



g(x)= ( --2)F(0,a;) + e 
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a primitive function of f{x) (see Appendix). When the 
slab thicknesses become infinitely large we obtain for the 
interaction between two bulk plates: 
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C. Sphere-slab geometry 



We next consider a third experimentally relevant sit- 
uation where object A is (again) a slab of thickness ca 
and density , while object _B is a sphere of radius R 



(5) 
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and density . L is the surface-surface distance of the 
two objects, while £ will be the center-to-plate distance 
{C = L + R) (see Fig.©. 



III. LOCAL FIELD EFFECTS AND PWS 




FIG. 3. Sphere-slab geometry. 

We compute the interaction between the two objects 
by summing the atom-slab result given by Eq.([3]) for each 
atom of the sphere B which lies at a distance d = C + r 
from the slab, with r G [—R, R]: 



UrX-'AC,eA)^ f dr {7:{R'-r')n^)urs'{C + r,eA) 

J-R 

(7) 



After several integration by parts, it is possible to show 
that the energy reads: 



(^^^ 

X {2uR [h {2u{C + eA + R)) + h {2u{C + ca ~ R)) 
- h {2u{C + R)) - h {2u{C - R))] 
~ [i {2u{C + eA + R)) - i i2u{C + ca - R)) 

-i {2u{C + R)) + i i2u{C - R))]} (8) 
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two successive primitives oi g{x) (see Appendix). In the 
limiting case of e^i — )■ 00, we obtain: 
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In Sec. Ill Al we calculated the PWS estimate of the 
atom-slab interaction by direct summation of the van 
der Waals atom-atom interactions. We now show that 
the same result can be derived within the scattering ap- 
proach, which will make it easier to compare analytically 
the PWS result and the exact Casimir interaction in the 
atom-slab geometry. 



Equivalence between PWS and the summation 
of reflection matrices 



The Casimir interaction at zero temperature between 
any two objects A and B at a distance R along z from one 
another can be expressed as follows within the scattering 
approach to Casimir forcesi^"— : 



UiR) 
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In this formula, TZa and TZ b are the reflection matrices of 
the two interacting objects, while the two e~'^^ factors 
arc propagation matrices. The term TZb e~'^^TZA e~^^ 
therefore corresponds to a complete roundtrip of the vac- 
uum field between the two objects A and B: the field 
propagates from B to A before being reflected by A, then 
it propagates back to B and is reflected by B. The inte- 
gration is on imaginary frequencies a; = icu of the electro- 
magnetic field, while the trace Tr runs over all indepen- 
dent field modes at a given frequency w = icu. Therefore, 
this formula takes into account the contribution of every 
electromagnetic mode to the Casimir effect. 

If one of the objects is an atom, only a small fraction of 
the field is scattered back to object B and we may safely 
neglect multiple reflections. Thus, if A or B - or both ~ 
are atoms, the general formula ([T0|) can be used at first 
order in TZb c^'^^TZa e"'^^, reducing to: 



U{R) = - 



he 
2^ 



du Tr hlBe-^^ nAe-^") 



(11) 



In particular, inserting the expressions of atom reflec- 
tion matrices (|13p into the general Casimir interaction 
formula at first order (|lip gives the van der Waals inter- 
action between two atoms (HJ. 

Let us now consider, as in Sec. Ill A[ the case of an 
isolated atom B in front of a slab constituted of atoms 
A (see Fig. [T|). Since formula (fTTI) is linear in TZa: it is 
equivalent to carry out a pairwisc summation of atom- 
atom interactions as in Sec. Ill Al or to sum the reflection 
matrices of the constituting atoms in order to obtain the 
reflection matrix of the slab. More precisely, if we write 
the slab reflection matrix as 7^., = 72,yi(i?^), where 
Ra is the position of each atom A of the slab (the isolated 
atom B being at the origin of coordinates), the atom-slab 
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interaction can be written using 



Ua-s{L) 



(12) 



We thus obtain the predicted result. It can also be un- 
derstood intuitively: summing the reflection matrices 
of the constituting atoms amounts to considering each 
atom of the slab as an independent scatterer. Each of 
these atoms then interacts independently with the iso- 
lated atom through the electromagnetic field in an equiv- 
alent manner to PWS. 



Calculation of slab reflection coefficients by 
summation of reflection matrices 



Having established that PWS is equivalent to summing 
the atom reflection matrices in the atom-slab geometry, 
we will now compare the exact reflection coefficients of 
a slab to the reflection coefficients obtained by summing 
the reflection matrices of the constituting atoms. For 
this, let us first calculate explicitly the latter coefficients. 

First, we write the position vector of an atom A of 
the slab as Ra ~ {rA,ZA = L + (a)- Thus (a ranges 
from to eA (see Fig. [1]), while ta € R^. In order to 
compute TZs = '^jyTi-A(rA,L + Ca), we calculate a ma- 
trix element of TZs in the plane wave basis. At a given 
frequency, each plane wave is characterized by its two- 
dimensional transverse (i.e. orthogonal to z) wave vector 
k, its longitudinal direction and its polarization p, trans- 
verse electric (TE) or magnetic (TM). This description 
is complete since the modulus of the longitudinal com- 
ponent of the wave vector can be deduced from the 

relation w = c\ 



\k^ 



Let us choose two plane waves (fci,pi) and (^2,^2), the 
first one with longitudinal component in the -I- direction 
along z and the second one with longitudinal component 
in the — direction along z. Their respective unit polar- 
ization vectors are therefore noted £^^(^1) and £~^{k2). 
The matrix element taken between our two plane waves 
of the reflection matrix of an atom A at (0, L) is^: 



{k2,P2\TiA\kl,Pl) 



v?aA{icu) 



2^260 



^P2 



(13) 



In this formula, K2 is defined as: |/c|| = ik2- 
We can write {k2,P2_\'TiA{rA, L + CA)\ki,Pi) = 
(fc2,P2|'7^A|fci,Pi)e"'(*^""^i)'''^e"(''^"'"i)f-*, taking into 
account the field's propagation between (0, L) and each 
atom A situated at (rvi, za = L + C,a)- 



We are now able to compute the matrix element of 
TZs = ^A '^Ai'TA, L + Ca) between the two plane waves: 



{k2,P2\TZs\ki,pi) = 2^(fc2,P2| TZAirA,ZA)\ki,Pi) 

A 

= n,{k2,P2\TZA\ki,pi) I drXe-^^--^'^-'^^ 



dC,Ae 



-(k2-Ki)Ca 



(14) 



Once the integrals are calculated and {k2,P2\TZA\ki,pi) 
is replaced by its explicit expression (jl3[) . it yields: 



{k2,P2\TZs\ki,pi) 



i2nfdik2-ki) (e-2-i-A _i)(i5) 

The scalar products of polarization unit vectors can be 
easily calculated by writing down the explicit form of 
these vectors as a function of k and k^—. Thus, e^, (fci) • 
epi(^i) = if pi 7^ P2, while eTsiki) ■ ^Tsiki) = 1 and 
^TAiiki) ■ ^TAiiki) = ^(1 + 2;^). We finally obtain: 

( k2 ,p2\TZs\ki,pi) 
= (-I)-(PO {2TT)H{k2 - h)Sp,p, rP'{h,u) (16) 
where o{TE) ~ 1 and (t{TM) = —1, and where: 

7rM^aA(zCu)^^_2Kie^ T 



" (4^6o)«? 

TE 



rl^^{k^,u) = r'^{k^,u)- [1 + 2--^ 



k^ 



(17) 
(18) 



Equation pB]) is totally general since it reflects the 
invariance laws of specular reflection: the transverse wave 
vector is conserved by specular reflection (this arises from 
the X and y invariance of the system through Noether's 
theorem), and TM and TE polarization are cigenmodcs 
of specular reflection. 

On the contrary, the reflection coefficients (fT7|) and 
psp are specific to our summation approach. Using 
these coefficients instead of the exact ones to compute 
the atom-slab Casimir interaction is equivalent to using 
the PWS approximation. We have proven this formally in 
Sec. nil A[ and it is straightforward to verify it explicitly 
by calculating the atom-slab Casimir interaction using 
formula pi|) with the form of TZs we have just found: 
the result obtained is exactly the one obtained by PWS, 
i.e. dSl). 



C. Comparison of the exact Casimir interaction 
and the PWS estimate in the atom-thin slab 
geometry 

In order to gain some insight into the error made by 
PWS, we now compare the refiection coefficients calcu- 
lated by summation in the previous section to the exact 
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ones. This is equivalent to comparing the exact Casimir 
interaction and the PWS estimate in the atom-slab ge- 
ometry. 

So as to obtain its exact reflection coefficients, the slab 
is considered as a quantum optical network obtained by 
piling up a vacuum/matter interface, propagation over a 
length CA in matter and a matter/vacuum interface. Its 
transfer matrix is therefore the product of the transfer 
matrices of the three elementary networks constituting 
it. For a given polarization p=TE,TM the slab reflection 
coefficients can be obtained from its transfer matris2i. 



sh 1] 



' slab 



sh(77 + ep) 



(19) 



In this formula, rj = e^iKm corresponds to propagation in 
the slab, with k„ = Vew^-^F, while 9'^^ = In (f^) 
and 6*^*^ = In (fj^f^) are linked to the Fresnel vac- 
uum/matter reflection coefficients through the relation 
rP ~ —e^^^. Here we have noted in = \k^\ the modulus 
of the longitudinal wave vector and /c = |fc| the modulus 
of the transverse wave vector. 

If the exact reflection coefficients ([TO]) are used to 
calculate the atom-slab interaction from formula (|lip . 
the expression of the exact Casimir-Polder interaction is 
obtained^S. 

As both exact ^ and PWS (jlTHlSp reflection coeffi- 
cients are nonlinear function of the slab thickness e^i, we 
compare them in the case of a thin slab (e^ — >■ 0). This 
limiting case is interesting since one reason generally in- 
voked to explain the errors of PWS is that it fails to take 
into account the screening of the electromagnetic field. 
If screening was the only problem of PWS, this method 
would be most accurate for a thin slab. 

A first order development in ca of the reflection coef- 
ficients obtained by summation (jlTHlSp gives: 



slab 



' slab V 



-~2tT Uy Ca 

2n Uy ca 



aA{icu) 
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1 + 2^ 



(20) 
(21) 



Similarly, the exact reflection coefficients (|19p become, at 
ffist order in e: 



^slab(^' ^) 



6^(6- l)u2 



2k 



2k 



(22) 
(23) 



We notice that, at first order in ba, the two TE re- 
flection coefficients pO)) and (|22p are identical only if 
e = 1 -f ATTTiy , while the two TM reflection co- 

efficients pT|) and ([23]) are identical only if = 2, that 
is to say if e — 1 vanishes. Thus, we can say that PWS 
gives exact results in the atom-thin slab geometry only 
if e = 1 + Anny with AnUy < 1, that is 



to say only in the case when the material of the slab is 
infinitely diluted. 



This condition found in the limit of thin slabs is very 
general. In fact it can proven to be sufficient for any slab 
thickness: the condition e(icu) = 1 + 47rn„ "i'''"^ with 

47r Hy "4^'^"'* ^ 1 is the first-order approximation of the 
Clausius-Mossotti relation: 



e{icu) — 1 An a{icu) 
e{icu) + 2 3 " 47reo 



(24) 



when Uy a{icu) is very small, i.e. at the diluted limit. If 
we insert this first-order approximation to the Clausius- 
Mossotti relation in the exact reflection coefficients (|19p . 
we obtain the reflection coefficients found by summation 
(|17lll8p for any slab thickness. 

Replacing the Clausius-Mossotti relation by its first- 
order approximation amounts to neglecting local field ef- 
fects. Therefore our analysis indicates that local field 
effects are crucial in the errors made by PWS. More- 
over, the fact that the condition of validity of PWS is 
not more restrictive for a thick slab than for a very thin 
one suggests that screening is not the most important of 
the effects that cause the errors of PWS. 



IV. INFLUENCE OF THE MATERIAL IN THE 
LONG-RANGE LIMIT 

The formal comparison of the Casimir interaction and 
its PWS estimate in the atom-slab geometry has enabled 
us to gain some insight into the reasons why PWS is not 
exact, emphasizing the role of local field effects among 
the different many-body effects that are not taken into 
account in PWS. However, such formal comparisons are 
quite difficult to carry out in general, and it is of great 
practical importance to know the magnitude of the error 
made by approximations such as PWS. 

The results that are generally cited about the errors 
of PWS involve perfect reflector o^^'^^ . In the long-range 
(retarded) limit, the ratio between the pairwise approxi- 
mation and the exact value of the atom-slab interaction 
is 1.15 if the slab is a perfect reflector. In the slab-slab 
geometry, this ratio is O.SOii. As the perfect reflector is 
an idealization, we may wonder how this ratio changes 
for real materials. In this section, we estimate numeri- 
cally the ratio between the pairwise approximation and 
the exact value of the Casimir interaction for dielectric 
materials in the atom-slab and slab-slab geometries. In 
order to achieve model-independent results we restrict 
ourselves to the long-range limit. 



A. Atom-slab interaction 

Let us in a flrst step suppose that the slab thickness ca 
is much larger than L. The influence of the slab thickness 
will be addressed later-on. 

Besides we suppose that the distance L between the 
atom and the slab is much larger than any relevant in- 
trinsic characteristic distance, such as the wavelengths 
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corresponding to the atomic transitions of the isolated 
atom, and any wavelength characteristic of the permit- 
tivity of the slab material. Practically, for a typical atom 
and a silicon slab with permittivity modeled by a Drude- 
Lorentz formula3^, this means that L > 1 /im. 

The exact atom-bulk plate interaction can be ex- 
pressed thanks to formula (fTTI) . where 72.p is the slab re- 
flection matrix (|16p with the exact reflection coefficients 
([T9)) . The latter reduce to the bulk Fresnel reflection co- 
efficients, since we are dealing with infinite slabs. The 
result reads2£: 



he 

TE 



2—-l\r 



TM 



(25) 



where r'^^ 



and r'^^' 



with K„ 



^y v?e{icu) -\- k^. In the long-distance limit, the propa- 
gation factor e~^'^^ in the atom-plate Casimir interac- 
tion formula psp vanishes except at small frequencies 
uj = icu, so that the polarizability of the isolated atom 
asiicu) and the permittivity e{icu) can be replaced by 
their static values Q!_b(0) and e(0). 

In the PWS approximation, the interaction will thus be 
given by formula (jj]). As before, the propagation factor 
g-2uL a^jiows for replacement of aA{icu) and asi^cu) by 
0^(0) and asiO) in this formula in the long distance 
limit. Explicit integration then shows that formula (|^ 
reduces to: 



rrPWS 

^ n — n 



23 hcn:^ aA{0)aB{0) 



(47reo)2 



(26) 



This long-distance limit of formula Q can be found 
by simple pairwise summation of the retarded Casimir- 
Polder interactionii. 

We notice that in the long-distance limit, q;b(0) is a 
simple multiplicative factor in both the exact Casimir in- 
teraction and its PWS estimate. Therefore, in this limit, 
the ratio U^^^ /Ua-p does not depend on any property 
of the isolated atom. Besides, in order to compare the 
PWS estimate to the exact Casimir interaction, Q!a(0) is 
expressed as a function of e(0) thanks to the Clausius- 
Mossotti relation 

It is now possible to compute the ratio U^^^ /Ua-p in 
the long range limit as a function of the only remaining 
material parameter e(0). The result obtained is traced in 
Fig.il 

We observe on Fig. [4] that the PWS method becomes 
correct when e(0) 1, that is to say when the bulk 
matter becomes infinitely diluted. This is consistent 
with the conclusions of Sec. IIIII Besides, we find the 
well-known result U^^^ /Ua-p = ~ 1.15 correspond- 
ing to the perfect reflector in the limit e(0) — >■ ocM^. 
An unexpected and interesting result is that the ratio 
U^^^ /Ua~p does not evolve monotonically between these 
two limits, but has a maximum at e(0) ~ 14.9, where 
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FIG. 4. Ratio of the PWS estimate of the atom-plate Casimir 
interaction and its exact value in the long-range limit. This 
ratio is traced as a function of the static permittivity. 



Ua^p^/Ua-p ^ 1.321. Moreover, this maximum cor- 
responds to realistic dielectric materials. For instance, 
we show on Fig. [4] that silicon is very close to it, with 
e(0) = 11.87 and U^^J^/Ua-p ^ 1.319. PWS thus makes 
a larger error for dielectrics such as silicon than for a 
perfect reflector if the atom-bulk plate interaction is con- 
sidered. The relative error can reach about 30%. 

We now consider a finite thickness ca for the slab and 
study its effect on the accuracy of the PWS method. We 
start from the derived expression ([3]) for the energy be- 
tween the slab and the atom, and then apply the long- 
distance limit so that the polarizability of the isolated 
atoms can be replaced by their static values. It leads 
to a result similar to the bulk case ([^S)) . except for a 
thickness-dependent factor: 



U. 



PWS 



(L,e^) = C/r-nL) 



L , 



(27) 



For the computation of the exact Casimir energy, we 
use Eq.(|25|) with the thickness-dependent Fresnel coef- 
ficients introduced in Eq. p9|l . The obtained ratio is 
traced in Fig. [5] for various values of the relative thick- 
ness e = ca/L which is dimensionless. We observe that 
the ratio of the PWS estimation over the exact result 
is also strongly dependent on the thickness of the slab 
compared to the distance between the two objects: while 
for a large thickness compared to the distance (e ^ 1) 
the bulk plate case is recovered for any static permittiv- 
ity, the ratio decreases dramatically when the thickness 
decreases. At the limit of infinite permittivity, or per- 
fectly reflecting objects, the ratio goes as expected to the 

thickness-dependent ratio ^1 — (1 -j- e^/i)"^^ ■ 
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FIG. 5. Ratio of the PWS estimate of the atom-slab Casimir 
interaction and its exact value in the long-range limit, for 
various values of the relative thickness e = ca/L. The case of 
a bulk plane is recalled with a dashed-line. 



B. Two slabs 

Let us now carry out the same study for two bulk ma- 
terials. In this geometry, the whole formula (fTU]) has to 
be used since multiple reflections arc not negligible. Us- 
ing equation (jl6[) for the reflection matrix of each plate, 
this formula gives the energy per unit area: 



he f°° 

Up-p{L) = — / du Tr In (l - 



TA e 



-2kL\ 



(28) 



This formula is only valid for specular reflection. The 
general formula (|10p for the Casimir interaction between 
any two objects was in fact first obtained as a generaliza- 
tion of the specific formula ([28]) ^. After writing explic- 
itly the trace over field modes, formula ([28|) becomes: 



TE\2^-2kL\ 



Up.p{L) = -4 / K [ In (1 - (r^^) 

47r Jo Ju 

+ ln(l-(r™)2e-2-^) ] (29) 



with the bulk rcfiection coefficients r and r such 
as defined after eqn. ([^5]) . Wc have assumed here that 
both plates are identical. In the case of perfect reflectors 
(r^ = 1), formula ([29]) reduces after some calculations to 

the original Casimir formula.^ Up^p{L) = —^f^j- 

In the PWS approximation, the interaction between 
two plates at any distance has been computed in Sec. |IT1 
formula ([6|). In the long distance limit the polarizabilities 
can be replaced by their static values, so that the general 
PWS result ^ can be simplified: 



23_ hcn^n^ aA{0)aB{0) 
120 



L3 



(47reo)2 



(30) 



This long-distance limit of formula ([6]) can also be found 
by simple pairwise summation of the retarded Casimir- 
Polder interactionii. Here we have a a = 0-3 a-nd n-^ ~ 



In order to compare the exact Casimir interaction 
and its PWS estimate, we use once more the Clausius- 
Mossotti relation ([M]) to express a{0) as a function of 
e(0). We now compute the ratio Up^^ /Up-p in the long 
range limit using the same method as in the previous sec- 
tion. The result obtained is traced as a function of e(0) 
in Fig. [SI 




since both plates are supposed to be identical. 



FIG. 6. Ratio of the PWS estimate of the plate-plate Casimir 
interaction and its exact value in the long-range limit. This 
ratio is traced as a function of the material's static permit- 
tivity. 

While the global aspect of the ratio UF^^ /Up-p in 
Fig. [S] is similar to the one plotted in Fig. g] (atom-bulk 
plate interaction), the error of PWS in the two plates 
geometry is found to be twice as large, reaching now up 
to 60%. 

Moreover, a striking feature in Fig. [6] is that the ratio 
Up^^ /Up-p is larger than one for most dielectric mate- 
rials (e(0) < 100), while it is smaller than one for the 
perfect reflector {U^^J^/Up^p = ^ ^ 0.797). Thus, 
the perfect reflector result does not provide a correct ex- 
planation as to the reasons of failure of PWS in most 
materials. This observation illustrates the complexity of 
the error made by PWS: the PWS estimate can be ei- 
ther smaller or larger than the exact Casimir interaction 
depending on the material for a given geometry. Further- 
more, screening of electromagnetic fields by slab atoms 
turns out not be the main responsible for the error as 
we would then expect the ratio Up^^ /Up^p to be larger 
than one for all values of e(0), which is not the case here. 
As we have argued in Sec. IIIH local field effects seem to 
play a crucial part in this error. 

We now consider finite thicknesses cajCb for the slabs 
and study the effect of these parameters on the accuracy 
of the PWS method. We start from the derived expres- 
sion ([5]) for the energy between the slabs, and then apply 
the long-distance limit so that the polarizability of the 
isolated atoms can be replaced by their static values. It 
leads to a result similar to the bulk case ([5^ . except for 



8 



a thickness-dependent faetor: 



X 1 



Us_s {L,eA,eB) = Up_p [L) 
1 



(1 + ^)^ 



[1 + ^ + ^] 



(31) 



For the computation of the exact Casimir energy, we 
use Eg. ip^ with the thickness-dependent Fresnel coefE- 
cients introduced in Eq. (|19p . For simphcity we will only 
consider the case where ba = gb- The obtained ratio is 
traced in Fig. [7] for various values of the relative thickness 
e ^ ca/ L = cb / L, as in the atom-slab case in Fig. [5l 
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FIG. 7. Ratio of the PWS estimate of the slab-slab Casimir 
interaction and its exact value in the long-range limit, for 
various values of the relative thickness e = ca/L = es/L. 
The case of two bulk planes is recalled with a dashed-line. 

We observe that the ratio of the PWS estimation 
over the exact result is also strongly dependent on the 
thickness of the slabs compared to the distance be- 
tween them: while for a large thickness compared to 
the distance (e ^ 1) the bulk plate case is recovered 
for any static permittivity, the ratio decreases dramat- 
ically when the thickness decreases. At the limit of 
infinite permittivity, or perfectly reflecting objects, the 
ratio goes as expected to the thickness-dependent ratio 

l-2(l + e)"V(l + 2e)" 



C. Sphere-bulk plate geometry 



between the sphere and the infinite slab L = £ — i? is 
much larger than any wavelength A characteristic of the 
permittivity e of the objects material. As in the previous 
cases the permittivity and polarizability in the reflection 
operators TZs and TZp can then be replaced by their static 
value e(0). 

In the PWS approximation, the interaction between a 
sphere and a plate has been computed in Sec. |TT1 formula 
([9]). In the long-range limit the general PWS result Q 
can be simplified to: 



C/, 



PWS 
sph- 



23 hc-nW' n^n^ 
30 {C^-R^Y 



aA(0)aB(0) 



(32) 



Here we have ua = O-b and 



since the materials 



for the sphere and the slab are supposed to be identical. 

In order to compare the exact Casimir interaction 
and its PWS estimate, we use once more the Clausius- 
Mossotti relation ([24]) to express a{0) as a function of 
e(0). We now compute the ratio Uf^^p/Ugph-p in the 
long range limit using the same method as in the previ- 
ous sections. The result obtained is traced as a function 
of e(0) in Fig. [HI for several values of the additional geo- 
metrical parameter 




FIG. 8. Ratio of the PWS estimate of the sphere-slab Casimir 
interaction and its exact value in the case of infinitely thick 
slabs in the long-range limit, for various geometrical parame- 
ters from 0.02 to 5000 (blue curves). This ratio is traced as 
a function of the static permittivity of the slab material. The 
black and red curves are a recall of the atom-slab and slab- 
slab cases, previously presented in Figs. |3]and|Sl respectively. 



We finish this study with the sphere-slab configuration, 
where a sphere of radius R has its center at a distance C 
from an infinitely thick slab. In this geometry the exact 
quantity for the Casimir energy Ugph-p has to be derived 
from the general scattering formula pO|) . with TZs and 
TZp the reflection operators on the sphere and on the 
plane, respectively. These operators can be expressed 
thanks to the bases of planar and spherical waves, as 
discussed in details in^^i^i. In the present study, we con- 
sider the long-range limit, in the sense that the distance 



The first observation is that the ratio U^^^p/Usph-p is 
dependent on the parameter : for small values (i.e. very 
large spheres) this ratio tends to the plate-plate result, 
recalled by the red curve. On the opposite, for large 
values of ^ (i.e. small spheres) the curves converge to the 
atom-plate result, presented with a dark curve, except 
for very large values of e. Indeed, the small-sphere limit 
{R <C L) for perfect mirrors yields a ratio || ~ 0.77, 
while the e(0) — >■ oo limit for an atom in front of a slab 
leads to a ratio ^ = 1.15. This shows that the two limits 
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{R ^ L) and (e(0) oo), of a small sphere and of a 
perfect reflector, do not commute. The ratio | is typical 
for this non-commutativitj«2i and can be traced back to 
the small parameter-limit for the Mie coefficient bg, or 
in other terms, to the impossibility to have a magnetic 
point-like dipole. 



providing excellent possibilities for discussions and ex- 
change. 

Appendix A: Expressions for the successive 
primitive functions 

The successive primitive functions involved in the 
derivations of the different energies are: 



V. CONCLUSION 

We have investigated the error made by the pairwise 
summation method in three geometries where the exact 
formula for the Casimir interaction is known, the atom- 
slab, the slab-slab and the sphere-slab geometry. The 
scattering approach to the Casimir effect has enabled us 
to show, through an analytical comparison of reflection 
coefficients, that the PWS result is incorrect even for an 
infinitely thin slab interacting with an atom, unless the 
material is infinitely dilute. This analysis has stressed the 
fundamental importance of local field effects among the 
many-body effects that PWS fails to take into account. 

We have then studied the influence of the material on 
the error made by PWS. This study has shown that the 
error made by the PWS method is much higher in the ex- 
perimentally relevant case of dielectric materials such as 
silicon than for perfect reflectors. We have reached this 
conclusion in the long range limit, both in the atom-slab 
and in the slab-slab geometry. The existence of a maxi- 
mum in the error made by PWS for the permittivity of 
usual dielectrics is not easy to understand intuitively and 
sheds light on the complexity of the error made by PWS. 
So does the fact that PWS can underestimate or overes- 
timate the Casimir interaction in the slab-slab geometry 
depending on the slab permittivity. Although the error 
made by PWS is influenced by the slab thickness, in both 
geometries studied, it turns out to be no more accurate 
in the case of thin slabs than in the case of thick ones. 
This result confirms that the fact that PWS does not 
take into account the screening of electromagnetic fields 
cannot be the only explanation of the errors made by this 
method. For silicon slabs, PWS was found to underesti- 
mate or overestimate the Casimir interaction depending 
on the thickness in both geometries, showing once again 
the complexity of the error made by PWS. 
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